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It is argued that the massive gauge field theory without the Higgs mechanism can well be set 
up on the gauge- invariance principle based on the viewpoint that a massive gauge field must be 
viewed as a constrained system and the Lorentz condition, as a constraint, must be introduced from 
the beginning and imposed on the Yang-Mills Lagrangian. The quantum theory for the massive 
gauge fieldis may perfectly be established by the quantization performed in the Hamiltonian or the 
Lagrangian path-integral formalism by means of the Lagrange undetermined multiplier method and 
shows good renormalizability and unitarity. 
It is the prevailing viewpoint that the massive gauge field theory can not be set up without introducing the Higgs 
mechanism. The first obstacle is the gauge-non-invariance of the mass term in the massive Yang-Mills Lagrangian for 
a massive gauge field. On the contrary, we present an argument to show that the conventional viewpoint is not true 
[1]. In fact, a certain massive gauge field theory can be well established on the basis of gauge-invariance principle 
£N| ■ without recourse to the Higgs mechanism. The basic ideas are stated in the following. 

(1) A massive gauge field must be viewed as a constrained system. In the previous attempt of building up the 
massive non-Abelian gauge field theory, the massive Yang-Mills Lagrangian density written below was chosen to be 
the starting point [2,3]. 
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where A® are the vector potentials for the non-Abelian massive gauge fields, 
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O : = d^Al - d v Al + gf ahc A\Al (2) 



are the field strengths and m is the mass of gauge bosons. The first term in the Lagrangian is the ordinary Yang-Mills 
Lagrangian which is gauge-invariant under a whole Lie group and used to determine the form of interactions among 
the gauge fields themselves. The second term in the Lagrangian is the mass term which is not gauge-invariant and 
only affects the kinematic property of the fields. The above Lagrangian itself was ever considered to give a complete 
description of the massive gauge field dynamics. This consideration is not correct because the Lagrangian is not only 
not gauge-invariant, but also contains redundant unphysical degrees of freedom. As one knows, a full vector potential 
Oh. A^(x) can be split into two Lorentz-covariant parts: the transverse vector potential At^(x) and the longitudinal vector 
potential A^{x), A^(x) — A^(x) + A^(x). The transverse vector potential A^(x) contains three independent spatial 
components which is sufficient to represent the polarization states of a massive vector boson. Whereas, the longitudinal 
vector potential A a ^ appears to be a redundant unphysical variable which must be constrained by introducing the 
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Lorentz condition 



& ; <p a = d»A° = (3) 

whose solution is A°^ = 0. With this solution, the massive Yang-Mills Lagrangian may be expressed in terms of the 
independent dynamical variables A^(x) 

C = -\F^F^ v + \m?A!FA^ (4) 

which gives a complete description of the massive gauge field dynamics. If we want to represent the dynamics in the 
whole space of the full vector potential as described by the massive Yang-Mills Lagrangian in Eq.(l), the massive 
gauge field must be treated as a constrained system. In this case, the Lorentz condition in Eq.(3), as a constraint, is 
necessary to be introduced from the onset and imposed on the Lagrangian in Eq.(l) so as to guarantee the redundant 
degrees of freedom to be eliminated from the Lagrangian. 

(2) The gauge-invariance should generally be required for the action written in the physical subspace. Usually, the 
gauge-invariance is required to the Lagrangian. From the dynamical viewpoint, as we know, the action is of more 
essential significance than the Lagrangian. This is why in Mechanics and Field Theory, to investigate the dynamical 
and symmetric properties of a system, one always starts from the action for the system. Similarly, when we examine 
the gauge-symmetric property of a field system, in more general, we should also see whether the action for the system 
is gauge-invariant or not. In particular, for a constrained system such as the massive gauge field, we should see 
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whether or not the action represented by the independent dynamical variables is gauge- invariant. This point of view 
is familiar to us in the mechanics for constrained systems. Certainly, in some special cases, the Lagrangian itself is 
gauge-invariant in the physical subspace so that the gauge-invariance of the action is ensured. This situation happens 
for the massless gauge fields and the massive Abelian gauge field. For the non-Abelian gauge fields, the infinitesimal 
gauge transformation usually is given by [3] 

SA a = D ab 0b (5) 

where 

D a J> = 5 a %-gf abc A^ (6) 

This gauge transformation is different from the Abelian one in that in the physical subspace defined by the Lorentz 
condition, i.e., spanned by the transverse vector potential A!j? , the fields still undergo nontrivial gauge transformations 

5A a T , = D%8 b (7) 

where 

D% = 5"% - gf abc A^ (8) 

Therefore, the mass term in the massive Yang-Mills Lagrangian written in Eq.(8) is not gauge-invariant. But, the 
action given by this Lagrangian is invariant with respect to the gauge transformation shown in Eqs.(7) and (8). In 
fact, noticing the identity: f abc A^f ' A b T)l = and the transversity condition (an identity): d^A^^ = 0, it is easy to see 

SS = -m 2 J d A x8 a d^A^^ = (9) 

This shows that the dynamics of massive non-Abelian gauge fields is gauge- invariant. Alternatively, the gauge- 
invariance may also be seen from the action given by the Lagrangian in Eq.(l) which is constrained by the Lorentz 
condition in Eq.(3). Under the gauge transformation written in Eqs.(5) and (6), noticing the identity f abc A afl A b l = 
and the Lorentz condition, it can be found that 



SS= -to 



J d i x9 a d^Al = (fO) 



This suggests that the massive non-Abelian gauge field theory may also be set up on the basis of gauge-invariance 
principle. 

(3) Only infinitesimal gauge transformations need to be considered in the physical subspace. In examining the 
gauge invariance of the action for the massive non-Abelian gauge fields, we confine ourself to consider the infinitesimal 
gauge transformation only. The reason for this arises from the fact that the Lorentz condition limits the gauge 
transformation only to take place in the vicinity of the unity of the gauge group. In other words, the residual gauge 
degrees of freedom existing in the physical subspace are characterized by the infinitesimal gauge transformations. 
This fact was pointed out in the pioneering article by Faddeev and Popov for the quantization of massless non- 
Abelian gauge fields [4]. Usually, the dynamics of massless gauge fields is described by the Yang-Mills Lagrangian. 
ft is well-known that the Yang-Mills Lagrangian itself is not quantizable, namely, it can not be used to construct 
a convergent generating functional of Green's functions even though it is gauge-invariant with respect to the whole 
gauge group. This is because the Yang-Mills Lagrangian contains redundant unphysical degrees of freedom and hence 
is not complete for describing the massless gauge field dynamics unless a suitable constraint such as the Lorentz 
condition is introduced to eliminate the unphysical degrees of freedom. In the article by Faddeev and Popov, the 
Lorentz condition is introduced through the following identity [4] 

A[A] J D(g)S[d^Af 1 ] = 1 (ff) 

which is inserted into the functional integral representing the vacuum-to vacuum transition amplitude. After doing 
this , the authors said that " We must know A [A] is only for the transverse fields and in this case all contributions 
to the last integral arc given in the neighborhood of the unity clement of the group" . The delta- functional in Eq.(f 6) 
implies d fl (A 9 )'^ = d^A® = which represents the gauge-invariance of the Lorentz condition because the Lorentz 
condition is required to hold for all the field variables including the ones before and after gauge transformations. In 
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the physical subspacc where only the transverse fields are allowed to appear, only infinitesimal gauge transformations 
around the unity element are permitted and needed to be considered in the course of Faddeev-Popov's quantization 
even though the Yang-Mills Lagrangian used is invariant under the whole gauge group. Obviously, there are no reasons 
of considering the gauge-transformation property of the fields in the region beyond the physical subspace because the 
fields do not exist in that region. By this point, it can be understood why in the ordinary quantum gauge field theories 
such as the standard model, the BRST-transformations are all taken to be infinitesimal. 

According to the general procedure, the Lorentz condition (3) may be incorporated into the Lagrangian (1) by 
the Lagrange undetermined multiplier method to give a generalized Lagrangian. In the first order formalism, this 
Lagrangian is written as [5] 

C = ^F a ^F* u - ^F aft " (dpA® - d v A1 + gf ahc A\Al) + ^m 2 A a ^A* + \ a d^A^ (12) 

where A^ and F* v are now treated as the mutually independent variables and A a are chosen to represent the Lagrange 
multipliers. Using the canonically conjugate variables defined by 

T\ a (r} — ^ _ pa , xa r _ t^kO = -^fc ' if ^/ = A; = 1 , 2 , 3 ; , . 

" ( ' ~ dA*» - + A ^° - { X a = -£ a , if /i - 0, {13 > 

the Lagrangian in Eq.(18) may be rewritten in the canonical form 

C = + A%C a - E$ip a - U (14) 

where 

C a = Q„ E a + g fabc A b E ck + ^ 



U = \(E a k f + \(F^f + \m 2 [(A a f + (Al) 2 ] (16) 

here E® = (Eft, E%) is a Lorentz vector, H is the Hamiltonian density in which F° are defined in Eq.(2). In the above, 
the four-dimensional and the spatial indices are respectively denoted by the Greek and Latin letters. From Eq.(14), 
it is clearly seen that the second and third terms are given respectively by incorporating the constraint condition 

C a = (17) 

where C° was represented in Eq.(15) and the Lorentz condition in Eq.(3) into the Lagrangian. 

Now, let us first perform the quantization of the massive non-Abelian gauge fields in the Hamiltonian path-integral 
formalism [5]. In accordance with the general procedure of the quantization, we should first write the generating 
functional of Green's functions via the independent canonical variables which are now chosen to be the transverse 
parts of the vectors A 1 ^ and E® 

Z[J] = ± J D(A a /, E a T »)ex V {i J d 4 x[E^A^ - H*(A a /, E a /) + J^A^}} (18) 

where Ti* (Ajf", E^l 1 ) is the Hamiltonian which is obtained from the Hamiltonian in Eq.(16) by replacing the constrained 
variables A a ^ and E^ 1 with the solutions of equations (3) and (17) 

H*(A a /, E a /) = U{A^, E a ») U =0 ,c»=o (19) 

As mentioned before, Eq.(3) leads to A a ^ — 0. Noticing this solution and the decomposition E atJ -(x) — E^ + E^(x), 
when setting 

=W(*) (21) 
where Q a (x) is a scalar function, one may get from Eq.(17) an equation obeyed by the scalar function Q a (x) 

K ab (x)Q b (x) = W a (x) (22) 

where 
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K ab (x) = S ab a x - gf abc A c ^{x)d x (23) 

and 

W a (x) = gf abc E b T »{x)A c Tll {x) - m 2 Af{x) (24) 
With the aid of the Green's function G ab (x — y) (the ghost particle propagator) which satisfies the following equation 

K ac (x)G cb (x-y) = S ab 8 4 (x-y) (25) 
one may find the solution to the equation (22) as follows 

Q a (x) = J d 4 yG ab (x-y)W b (y) (26) 

From the expressions given in Eqs.(21) and (26), one can see that the E^(x) is a complicated functional of the 
variables and Ej? so that the Hamiltonian TL*{A C ^ 1 , E^) is of much more complicated functional structure which 
is not convenient for constructing the diagram technique in the perturbation theory. Therefore, it is better to express 
the generating functional in Eq.(18) in terms of the variables A® and E®. For this purpose, it is necessary to insert 
the following delta-functional into Eq.(18) [5] 

5[A a »]5[E a L » - E a L "(A^,E^)] = detM <5[C Q ]% Q ] (27) 

where M is the matrix whose elements are defined by the Poisson bracket 

M*{z,y) = {<?■(*), p»(y)} =fd*x{ s ie a (x) J?; {x) - j^j^} m 
= Df{x)d^{x-y) 

The relation in Eq.(27) is easily derived from equations (3) and (17) by applying the property of delta- functional. 
Upon inserting Eq.(27) into Eq.(18) and utilizing the Fourier representation of the delta-functional 

5[C a ] = J D( V a /2iryJ d4xr i aca (29) 

we have 

Z[J] = i / D(A«, V a )detMS[d^A-} exp{z / d 4 x[E^A* 

+rfC a - H{A a ^, E a ^) + J a M£]} y ' 

In the above exponential, there is a i?Q-related term Eq(8oAq — dor] a ) which permits us to perform the integration 
over Eq, giving a delta- functional 

S[d A a - d v a ] = det\dv\-H[A a - v a ] (31) 

The determinant de^c^l -1 , as a constant, may be put in the normalization constant TV and the delta- functional 
5[Aq — rj a ] will disappear when the integration over if is carried out. The integral over E% is of Gaussian- type and 
hence easily calculated. After these manipulations, we arrive at 

Z[J] = j ! jD(A-)detMd[d^A-}exp{ifd i x[-\F^F- I/ 

+ i m 2 i4 «M° + J°MJ]} [ ' 

When employing the familiar expression [4] 

detM = J D (c\C a )e f dixdivCaix)Mabix - v)cb(y) (33) 

where C a (x) and C a {x) are the mutually conjugate ghost field variables and the following limit for the Fresnel 
functional 

5[d»A a u ] = lim C[a]e^/ ^-A^ (34) 



4 



where C[a] ~ Ila^;^) 1 ^ 2 anc ^ supplementing the external source terms for the ghost fields, the generating functional 
in Eq.(32) is finally given in the form 

Z[J,l£] = ± J D{A^C a ,C a )exp{i J d i x[C eSf + J^Al+Tc a + C a i a ]} (35) 

where 

C eff = --F a » v F« v + \m 2 A a ^A a il - -L(o>^) 2 - d^C a D ab C b (36) 

which is the effective Lagrangian for the quantized massive non-Abelian gauge field in which the third and fourth 
terms are the so-called gauge-fixing term and the ghost term respectively. In Eq.(36), the limit a — > is implied. 
Certainly, the theory may be given in arbitrary gauges (a ^ 0). In this case, as will be seen soon later, the ghost 
particle will acquire a spurious mass /i = y/am. 

To confirm the result of the quantization given above, let us turn to quantize the massive non-Abelian gauge fields 
in the Lagrangian path-integral formalism. For later convenience, the massive Yang-Mills Lagrangian in Eq.(l) and 
the Lorentz constraint condition in Eq.(3) are respectively generalized to the following forms 

C x = --F^F« V + ±m 2 A^A« -\a{\ a ) 2 (37) 

and 



d»Al + a\ a = (38) 

where X a (x) are the extra functions which will be identified with the Lagrange multipliers and a is an arbitrary 
constant playing the role of gauge parameter. According to the general procedure for constrained systems, Eq.(38) 
may be incorporated into Eq.(37) by the Lagrange undetermined multiplier method, giving a generalized Lagrangian 

C x = -\F a ^F^ v + ^m 2 A a ^Al + X a d^A a ^ + ^a(\ a ) 2 (39) 

This Lagrangian is obviously not gauge- invariant. However, for building up a correct gauge field theory, it is necessary 
to require the action given by the Lagrangian in Eq.(39) to be invariant with respect to the gauge transformations 
shown in Eqs.(5) and (6) so as to guarantee the dynamics of the gauge field to be gauge- invariant. By this requirement, 
noticing the identity f abc A afl A b 1 = and applying the constraint condition in Eq.(38), we have 



SS X = / (fxd v A a Jx)d lx {Vf{x)6 b {x)) = (40) 

a ' 



where 



Vf(x)=6 ab £-d* + Df(x) (41) 

in which /i 2 = am 2 and D ab (x) was defined in Eq.(6). From Eq.(38), we see \d v A a v = -A a ^ 0. Therefore, to ensure 
the action to be gauge-invariant, the following constraint condition on the gauge group is necessary to be required 

d%{Vf{x)9 h {x)) = (42) 

These are the coupled equations satisfied by the parametric functions 9 a (x) of the gauge group. Since the Jacobian 
of the equations is not singular, detM ^ where M is the matrix whose elements are 

ivi yx,y) - seb{y) I0=o 

= S ab (n x + ^ 2 )S%x - y) - gf abc dg(Al(x)6 4 (x - y)) 

the above equations are solvable and would give a set of solutions which express the functions 9 a (x) as functional of 
the vector potentials A^(x). The constraint conditions in Eq.(42) may also be incorporated into the Lagrangian in 
Eq.(39) by the Lagrange undetermined multiplier method. In doing this, it is convenient, as usually done, to introduce 
ghost field variables C a (x) in such a fashion [3-5]: 9 a (x) = <;C a (x) where <r is an infinitesimal Grassmann's number. 
In accordance with this relation, the constraint condition in Eq.(42) can be rewritten as 
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dHP£C b ) = (44) 

which usually is called ghost equation. When this constraint condition is incorporated into the Lagrangian in Eq. (39) 
by the Lagrange multiplier method, we obtain a more generalized Lagrangian as follows 

C x = -jF a ^F^ + ^m 2 A a>i A^ + \ a d^A^ + ^a(\ a ) 2 + C a d»(VfC b ) (45) 

where C a (x), acting as Lagrange multipliers, are the new scalar variables conjugate to the ghost variables C a (x). 

At present, we are ready to formulate the quantization of the massive gauge field in the Lagrangian path-integral 
formalism. As we learn from the Lagrange multiplier method, the dynamical and constrained variables as well as the 
Lagrange multipliers in Eq.(45) can all be treated as free ones, varying arbitrarily. Therefore, we are allowed to use 
this kind of Lagrangian to construct the generating functional of Green's functions 

Z[J">,?,?] - ^jD(A-,C a ,C a ,\ a )exp{iJd 4 x[C x (x) 

+J a K(x)Af l (x)+£, a (x)C a (x) + C a (x)£ a (x)}} [ ' 

where D(A®, ■ ■ ■ , X a ) denotes the functional integration measure, J", £° and £ a are the external sources coupled to the 
gauge and ghost fields and N is the normalization constant. Looking at the expression of the Lagrangian in Eq.(46), 
it is seen that the integral over X a (x) is of Gaussian- type. Upon completing the calculation of this integral, we finally 
arrive at 

Z[J<w,T,?] = ^SD(A^C a ,C a ,)exp{iJd l x[C eff {x) 

+J a »(x)A«(x)+£, a (x)C a (x) + C a (x)e(x)}} [ ' 

where 

Ceff = -]F a ^F^ + \m 2 A a »A a - ^-(^A") 2 - d< x C a VfC b (48) 

is the effective Lagrangian given in the general gauges. In the Landau gauge (a — > 0), The Lagrangian in Eq.(48) just 
gives the result in Eq.(36). It has been proved [1] that the above quantization carried out by means of the Lagrange 
multiplier method is equivalent to the Faddeev-Popov approach of quantization [4] . 

There are three points we would like to emphasize: (1) In the quantization by the Lagrange multiplier method, 
the gauge-invariance is always to be required even in the arbitrary gauge. Moreover, it has been found that the 
action given by the Lagrangian in Eq. (48) is invariant under a kind of BRST-transformations [6] . Thus, the quantum 
non-Abelian gauge field theory is set up from beginning to end on the firm basis of gauge-invariance. (2) In the 
Lagrangian path-integral formalism, as shown before, the quantized result is derived by utilizing the infinitesimal 
gauge transformations. This result is identical to that obtained by the quantization in the Hamiltonian path-integral 
formalism. In the latter quantization, we only need to calculate the classical Poisson brackets without concerning 
any gauge transformation. This fact reveals that to get a correct quantum theory in the Lagrangian path-integral 
formalism, the infinitesimal gauge transformations are only necessary to be taken into account and thereby confirms 
the fact that in the physical subspace restricted by the Lorentz condition, only the infinitesimal gauge transformations 
are possible to exist. (3) From the generating functional shown in Eqs.(47) and (48), one may derive the gauge boson 
propagator as follows 

which is of good renormalizable behavior. In the zero-mass limit, this propagator with the massive Yang-Mills 
Lagrangian and the generating functional together all go over to the results given in the massless gauge field theory, 
different from the quantum theory established previously from the massive Yang-Mills Lagrangian alone without any 
constraint [7-10]. For the previous theory, there occurs a severe contradiction in the zero-mass limit that the massive 
Yang-Mills Lagrangian in Eq.(l) is converted to the massless one, but, the propagator is not and of a singular behavior. 
In particular, the previous theory was shown to be nonrcnormalizable [3, 7-9] because the unphysical longitudinal 
fields and residual gauge degrees of freedom are not excluded from the theory. 

Up to the present, we limit ourself to discuss the gauge fields themselves without concerning fermion fields. For 
the gauge fields, in order to guarantee the mass term in the action to be gauge-invariant, the masses of all gauge 
bosons are taken to be the same. If fcrmions are included, Obviously, the QCD with massive gluons fulfils this 
requirement because all the gluons can be considered to have the same mass. Such a QCD, as has been proved, is 
not only renormalizable, but also unitary [6] . The renormalizability and unitarity are warranted by the fact that the 
unphysical degrees of freedom in the theory have been removed by the constraint conditions in Eqs.(38) and (45). 
The gauge- fixing term and the ghost term in Eq.(49) just play the role of counteracting the unphysical degrees of 
freedom contained in the massive Yang-Mills Lagrangian as verived by the perturbative calculations [6]. 
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